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Abstract. We develop a rigorously controlled multi-time scale averaging 
technique; the averaging is done on a finite time interval, properly chosen, 
and then, via iterations and normal form transformations, the time intervals 
are scaled to arbitrary order. Here, we consider as an example the problem of 
finite dimensional conservative dynamical system, which is quasiperiodic and 
dominated by slow frequencies, leading to small divisor problems in perturba- 
tive schemes. 



1. Introduction 

The aim of this work is to initiate the multiscale time-averaging analysis of 
repeated normal form finite time interval averages. That is we study the problem 
of general dynamical system type, written in the form 

(1.1) i^c = (3A(t)c 

where c is restricted to finite dimension, we choose A(t) to be a hermitian matrix, 
for each t, that is we restrict ourselves to conservative systems , i.e., 

N 

(1.2) ||c(t)|| 2 =^| Cl W| 2 = || C -*(0)|| 2 . 

i=l 

where N is the dimension of space. We construct finite time interval averages of 
A(t) 



(1.3) A { n) = - I A(s)ds 

J JnTo 

as approximation of the dynamics on a finite time interval depending on Tq. The 
equation is solved for the piecewise constant evolution, and we peel off the approx- 
imate solution from the exact, and derive an equation for the "leftover". 

Then, we introduce a normal from transformation, based on the method of [8j [9] ; 
this is an almost identity transformation of the "leftover". We then show, following 
(properly modified) arguments of [51 [3], that the resulting equation for the new 
quantity, satisfies the same equation as (|1.1|) . but with 



1 



(n+l)T 



3/2 



(1.4) p^p 



(1.5) A(t) -> A(t) 

with A(t) containing only "high frequency" terms. By repeating the above iteration 
we get, to any order, an operator V n (t), which solves to order /3 3 ™/ 2 . namely, the 

l 
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leftover, cm satisfies, 

(1.6) 

with 

(1.7) 



c(t) = V X V 2 



V M -lCM(t), 



V n — U\ n U2,nU n , 



where Ui iTl and V^n ar e unitary "peel off" transformations while U~ x is a nor- 
mal form transformation. The normal form transformation is an almost identity 
transformation, 

(1.8) Un = I + O (/3 3 "/ 2 ) . 

The cm (t) satisfy the equation similar to (|l.lll 

(1.9) 



8 



i^c M {t)=^ M/2 A M {t)c M {t). 

For insuring that the new A^> are hermitian after each iteration, we need to 
modify the normal form structure of [5J[S], and apply the following : 

Proposition 1. Let, for each t > the vector family c(t) € C N satisfy the follow- 
ing: 

(1.10) 11^)11 = 1, 



(1.11) 



c(t) = U(t)c\0), 



(1.12) 

with U linear. Furthermore 
(1.13) 

Then, 
(1.14) 
with 



U(tf = u(ty\ 



< oo. 



i^M=A(t)U(t), 



dt 



At = A. 



dc 



Proof. From 

and the conjugate equation, we get that 

d 

(1.16) 



i^rc(t) = A(t)c{t) 



i— c*c = c ■ A(t)c- A(t)c ■ c — 0. 



Hence , for all c(t), (c(t), A(t)c(t)) is real valued. Since U(t) is Unitary, this holds 
for ALL vectors in the space. Hence A(t) is hermitian for all t. 

The method of the present work involves an infinite hierarchy of averaging ap- 
proximations resulting in the effective reduction of f3 via (|1.4|l , allowing the validity 
of the approximation over increasingly longer time-scales, Tq. This is a substantial 
improvement compared to known methods where averaging in one stage is used as 
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described in [7] (see p. 379, App. E 2.2 and p. 390) and (4). The present work is an 
important step in this direction proposed in [7]. Finally, we demonstrate how the 
method works for the case of almost periodic A{t): 

oo 

(1.17) A(t) = J2 A ^ t , 
with 

(1.18) \\AjW < c(jy a , for some, a > 1 

and the interesting case wj — > 0, as j — > oo . Here Aj are constant N x N hermitian 
matrices. The interest in these kind of systems stems from the fact that many dy- 
namical systems can be modeled and or approximated by such equations. A related, 
but different example is the nonlinear system coming from Nonlinear Schrodinger 
Equation (NLSE) with a random potential term [3]. It leads naturally to a 
system with 



;i-i9) A n , n {t) = Y J v^ k i 



1,3 

where the regime 

(1.20) 0</3<l, 

is of great interest [T] . More generally, Hamiltonian dynamical systems with Hamil- 
tonians of the form 

(1.21) H = H Q + (3H 1 

can be started in this way. By solving exactly for the dynamics generated by Ho, 
we are reduced to a system similar to (jl.l[) . □ 

2. Averaging 

In this section an averaging of the matrix A (t) will be introduced. This averaging 
can be performed successively. Equation (|1.1[) is replaced to a hierarchical set of 
equations where A (t) and c(t) are replaced by A n and c n . The starting point of 
the hierarchy is Aq = A and cb = c satisfying the basic equation 

(2.1) i^c = 0A(t)c 

where c = (c^ ...c^) is a vector and A(t) is a matrix. We assume the following: 
Assumption 1. A(t) is, for each t, hermitian N x N matrix, satisfying, 

(2.2) sup 11^(4)11 < M < oo. 

t 

here, I denotes derivative with respect to t. /3 is a small positive number < /3 <C 1. 

Under Assumption [TJ the solution of (|l.ll) exits and is bounded uniformly in the 
usual vector norm on C™ in time, since, due to the hermitian property of A(t), 

(2.3) c{t)=Te- l0 K A{s)ds c{O) 

has a constant norm. T stands for the time ordering. Moreover, by taking deriva- 
tives, we get: 



(2.4) sup||?(t)||=sup 

i i 



-A{t)c{t) 

i 



< p svp \\A(t)c(t)\\ < pM 
t 
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We introduce now averaging over intervals of size Tq. Define the average matrix 



(2.5) 



I(n) _ j_ 
A Q — — 



1 



(n+l)T 



A(s)ds. 



To JnTa 

At a latter stage the relation between To and j3 will be defined. The global averaged 
function is 



(2.6) 



{n) 



for 



nT a <t< T (n + 1). 



Now one defines the propagator related to the averaged A, 



(2.7) U (t) = g-WCt-^o) . . . e -ipA 0) To 
and its inverse 



for 



nT Q < t < (n + 1)T 



(2.8) Uq 1 ^) = e +^4 0)T ° e +^o 1)r o . . . e ipA™(t-nT ) 
It is easily verified that: 

(2.9) U (t)U 1 (t) = l Uo 1 (t)U (t) = l. 
By direct differentiation one finds: 



(2.10) 



d 



i^-Uoit) = /3A^Uo(t) = /3A 9 (t)U (t), 



while the inverse satisfies 
d 



(2.11) 



Therefore, Uq is the propagator of the partially averaged motion generated by Aq (t) . 
The corresponding solution is, 

(2.12) £ (t) = UoHmt)- 

It satisfies 



(2.13) 



d _ 
l di Cl 



4t U ^ 



= [-pU^Alit)] c(t) + U^ 1 [PA{t)m] 
Where ([2~T0| and (|2~TT|) were used. Using (|2~12|) one gets 
d ^ 



^ci(t) - -pU^mt)U a {t)c x {t) +/3U 1 A(t)U (t)c 1 (t) 



d 



(2.14) 
Hence 

(2.15) i^ci(t) = PUo 1 [A(t) - A 9 (t)] tf (*)ci(t) 
This equation is analogous to (|1.1[) if written in the form 

(2.16) 
with 

(2.17) A 1 (t) = Uo 1 [A(t)-A 9 (t)]U (t) 



=PA 1 (t)3 1 {t) 

at 



Performing on (|2.16j) operations similar to the ones performed on (|2.1[) will lead to 
the next stage of the hierarchy. Before doing that we analyze the meaning of the 
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transformation to (|2. 16|) . The dynamics generated by A\ (t) is the peeling of that 
of A(t) by the average dynamics 



i An+1)T Q 

A s (t) = A^([t/T }) = ¥ I MsUls. 



JnT 



(2.18) - A(n} ' 
where, 

(2.19) n = [t/T ] , 

is the integer part of i/Tn. We turn to estimate the quantity 



{A{t')-Al{t'))dt' 
(2.20) 
Due to and (|23)l 
(2.21) 
Therefore, 



r/f' 



A(f)-4 0) (O 



2'ii, 



f/f' 



A(i')-4V) 



f/f' 



t-nTo 



A(t')-A^(t') 



(fe+l)T 



feT 



df {A(t>)-A°(t>))) = 



dt> (A(t)-A%(t)) 



< max t < 



A(t')-4 n V) (t-nTo) 



(2.22) 



< max 
t' 



A{t')-A^\t>) 



Tn 



where nTn < £' < (n + l)To. This result can be summarized as: 
Lemma 1. 



(2.23) 



{A{t')-Al{t'))dt' 



< 2 



tT 



Here \\\A\\\ t = sup < 4 ,< f and we used that |||^ n) (t)||| t < |||A(t)||| t . 

In what follows we choose 



(2.24) 



Tn 



1 

TP 



In order to estimate the variation of A\ it is useful to estimate the quantity: 



(2.25) 

or explicitly 
(2.26) 



h= I dt'A^t 1 ) 
Jo 



h= I UoHt')[A(t')-A 3 (t')]Uo(t')dt' 
h 
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For this purpose we note that integration by parts for matrices Mi is: 



Mi{t')M 2 (t')M 3 (t')dt' = Mi{t) 



M 2 {t')dt' 



M 3 (t) 



M 2 (s)ds 



M 3 (t') 



M 2 {s)ds 



dM 3 (t') 



dt> 



(2.27) - / dt'Mxtf 

Jo 

Therefore, using p3g]) . (j2~T0| and (j2TTj) . 

(2.28) h = U^{t)^[A{t')-Al{t')]dt^U Q {t) 
(PU^A* (f )) ^ [A( S ) - ig( s )] L/ (t'K 

Using (|2.23ll and ()2.24p one finds that the first term is of order, 
(2.29) 

while the other two terms are of order 



1 

VP 



(2.30) 
So that 
(2.31) 



In all these estimates we use the bound (|2.22[) that is independent of p. We turn 
now to the higher levels of the hierarchy. First, we introduce the global average of 
Ax (as we did for A = A ) and the transformation U\ corresponding to C/ - Let us 
define: 



X r(n+l)T 

A[ n) = I A x {t)dt 



nT 



(2.32) 

Where Ai(t) is given by (|2.17p . In analogy with (|2.6p we define 
(2.33) A{ (t) = A[ n) for nT Q <t<{n+ 1)T . 
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To estimate this quantity we apply the integration by parts (|2 .27|) to (|2.17l) . 



A™ = - 

1 rri 



(2.34) 







To J 


nTg 


1 

Th 


u \t) ( 






To J 


nTo 




r(n+l)T 


To J 


nT 



U - X (t) [A(t) - A{(t)]U (t)dt 
[A{t')-Al{t')}dAuo{t) 



(n+l)T 
t=nT 



t' 



dt 
Uo X (t') 



(A(s) - A 9 (s))ds 



U (t>) 



(A(s) - A°(s))dt 



-U (t')dt> 



The first term is zero since at time that is an integer multiple of T , /„ [A(t) ~ A 9 Q (t)] 
vanishes (see definition (|2.33[) ). From (|2 . 23[) and (|2.24[) we conclude that 



(2.35) 



A(s)-A g (s) 



< 



2kA 

VP' 



where is the norm of A. The derivatives dU^ 1 /dt and dUo(t)/dt are of order /3 
due to (|2.10p and (|2.f ip . therefore the second and third terms are of order y/]3. In 
conclusion, 

(2.36) \Al(t)\=o(y/p 



f7 _1 (f) = e l0A ^ )To ■ ■ ■ e ^A n) (t-nT ) 



Define in analogy to (f2~7]) and ([2~8]) . 

(2.37) C/i(i) =e- i ^ n) (*- nT °)---e- 

and 
(2.38) 

In analogy to (12.10p and (|2.1ip one finds 



(2.39) 

and 

(2.40) 



i^U x {t)=pA{{t)U x (t) 



In analogy to (|2.12p define now 

(2.41) c 2 (t) = U^ 1 (t)c 1 (t) 

and develop an equation analogous to (|2.16p . namely to 



(2.42) 



at 



For this purpose we follow the steps (|2.13p - (|2.1(jp . 



(2.43) ij^(t) 



c^ + U^it) 
pU^A{{t)U' 2 c 2 {t) + pU^Al (t)U^(t)c 2 (t) 



(2.44) 



i-c 2 (t) = pU-^t) [A^t) - A{{t)] U^c^t). 
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Define 
(2.45) 
leading to 

(2.46) 



A 2 (t) = Ur 1 [A 1 (t)-Al(t)]U 1 (t), 

z|-c 2 = (3A 2 (t)c 2 (t). 
ot 



Equation (|2.46| is similar in nature to (|2. 16|) and is the following equation in the 
hierarchy. Estimates of A 2 can be performed in the same way as the estimates of 
A\, The integral 

(2.47) I 2 = I dt'A 2 (t') 

Jo 

can be estimated in the same way as Ii in (|2.28l) . Namely, 

(2.48) h = U{\t) (j [Ai(f) - A{(t')] dA U (t) 

- j* {pU?A{ (t')) N* [Axis) - A{(s)] dA U^dt' 



Due to ([236) 
(2.49) 



fj ! !/' ) ( / [Axis) - A\ (a)] d'tj (-i/3lf (*')) Ui(t')dt'. 



A 



0) 



while by reasoning similar to (I2.22p 



(2.50) 

and 

(2.51) 



(A t (f) - A[ 9) (t')) dt' 



O 



h = [^)- <)t,n 



VP 



To estimate the magnitude of A 9 defined (in analogy to ()2.33|) ) as 



(2.52) 
where 

(2.53) 



A 



(a) ^(") 



AV = 



1 



for nT < t < (n + 1)T , 

(n+l)T 



Tn 



JnTo 



A 2 (*') dt' 



We now repeat the calculation of (|2.34[) . and use (|2.40[) and (|2.41|l combined with 
(|2.36p to estimate dU\/dt and dU^ 1 /dt. These derivatives are of order /3 3 / 2 . Con- 
sequently, 



(2.54) 



.4 



(</) 



0(/3). 



We have constructed explicitly the first two stages of the hierarchy: 

(2.55) c= U ci 
and 

(2.56) ci = Utc 2 . 
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The process can be continued further, defining 
(2.57) A n+1 (t) = U-\t) [A n (t) - A\{t)] U n (t). 

The equations similar to (|2.16[) and (|2.46[) are 

d 



i-^Cnit) = PA n (t)Cn (t). 



(2.58) 

The results of the present section can be summarized in 
Proposition 2. 



(1) 

(2) 
(3) 

(4) 
(5) 

(6) 



C= UqUi ■ ■ ■ U n -iCn 



A n {s) - AUs)ds 



< 



const 



\A n (s)-AUs)\=0(l) 
A\ = (>/ 2 ) 



In+l 



U-\ S ) [A n (s) - Ai(s)] U n (s) < O + O (>/ 2 i) 



(1) is a generalization of (|2.16[) and (|2.46p . 

(2) results of a repeated application of the transformation like (|2.55p and (|2.56p . 

(3) is a generalization of (|2.50p . 

(4) is a generalization of (j2"^T|) - (|2~2"3"|) . 

(5) is a generalization of (|2.36|) and (|2.54[) . 

(6) is a generalization of (|2.3ip and (|2.5ip . 

Remark 1. If the limiting operator exists, we get that is solves the following self- 
averaging equation: 



(2.59) 
where 



1 



To 



U-\s) [4>o0) - Q(s)] UMds = Q 



(2.60) 



Uoo = lim U n Aoo = lim A n Q = lim A 9 n . 
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3. Normal form transformation 

In this section a normal form transformation will be applied to c 2 (t) resulting 
in a new vector (t) . This vector satisfies an equation similar to (|f .ip but with (3 
replaced by /3 3 / 2 . We start from (|2.44|) 

(3.1) i^c 2 = pU{\t) [A x (t) - A{(t)] U x {t)c 2 {t)- 

Integrating by parts, we get: 

rt 



c 2 (t) - 3,(0) = -ij3 / U^it') [A, - Af(t')] U^ci^dt' 
Jo 

(3.2) = -ipiJ^it) U [Ai(0 - A{{t')} dt'^j U^Mt) 

+ if3 f* dt' (f* [At^-Afis^dsjm (t)c 2 (t) 

+ if3 j\t'U^{t>)(j^ [ Al {s)-Al{s)\d^^-{t')Ut') 
+ f dt'U^it') ( f [A^s) - A{{s)] ds) Uxit')-^^') 



10 \J0 

Using (HSU), ([2^0)1 and one finds 



(3.3) c 2 (t) - c 2 (0) + ipU^(t) (j [A^t') - A{{t')] dt'^j EM*)3i(t) 

= z/3 f dt 1 (+i(3U^(t)) A{{t) ( f [A x {s) - A?(s)] ] U x {1>)c 2 {t') 



+ i/3 / dt'f/fV) / ^(a) - Al{a)] ds (-i/3)lf (t')C/i^)c 2 (i') 



+ %p ^ dt'U^(t') ^ ^ ^(a) - Afia)] ds J U 1 {t l )(-iP)A 2 (t')c 2 (t l ) 
= -f3 2 f dt'U^{t')A{{t>) f [A^a) - Alia)} EMf)3»(f) 



•/(] 



+ /3 2 f\t'U^{t') [Axis) - Alia)] dsj AUt'W^it') 



ct / rt' 

p2 / J+lTT-lUI\ 







+ /3 2 ^ df'C/f 1 ^') y / - Af(s)] ds ) U l it')A 2 {t')c 2 it'). 

By (|2~50l) and Proposition 



(3.4) 



(Ai(«)-if(*)d. 



< O 



1 

VP 
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and bv 12361 

(3.5) \A 2 \~0(1) and A\ ~ 0(y/p), 

while A\ = (j3) by (|2.54p . The first two terms on the RHS are of order j3 2 and the 
last is term of order /3 3 / 2 . We turn now to perform the normal form transformation. 
For this purpose we rewrite (13.3P in the form 

[Atit') - A{{t')] dt'Uf(t)) Ui (t)c 2 (t)~c 2 (0) = RHS. 



(3.6) c 2 (t)+i/3U^(t) 

Then we define, 
(3.7) 



TfJ — 



Uc 2 



U = l + if3U^(t) / - A{(t')} dt' U^t) 



where 
(3.8) 

In these terms (|3.6j) takes the form 

(3.9) (t) = Uc 2 (t) - c 2 (0) = RHS, 

and differentiation of this equations with respect to time yields, 

. d 



(3.10) 
where 



dt 



A(t) 



+ UrHt) I [A^-Al^U^A^t) 



U^{t)A{(t)) (J [A^s) - A{(s)] ds ) U^t) 
[A^s) - A{(s)] ds) A{{t)U x {t) 

U~\t). 



(3.11) 



By (|3.4p and (|3.5p and A(t) = O (1). From its definition it is clear that U is also 

of order 0(1). It should be noted, however, that the normal form transformation 
above, is not unitary. Therefore, there is more than one way, by adding a correction 
to A, of higher order in /3, to make it a unitary operator. We do that by redefining 
U: 

Where B(t) = J„ [Ai(i') — A\ (t')\ dt' is self-adjoint. By choosing U to generate the 
normal form transformation, we insure that the L 2 norm does not change in time. 
Therefore, by Proposition [TJ the overall flow, generated by a product of unitary 
flows, is given in terms of a time dependent self-adjoint generator. 
Equation (|3.10p is similar to (jl.ip with the replacement, 



(3.12) 



c 
A 



A 

^3/2 _ 
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This reduces j3 effectively and increases the averaging time To. We can continue 
the process and find c n (t), such that 

(3.13) c=V 1 ---V n ^ 1 c n . 
Using (f2~55|) and ((2~56|) 

(3.14) c = f/iC/ 2 c 2 , 
and using ()3.7[) 



(3.15) c 2 = ET 1 *' 



2 ■ 



If one defines, 

(3.16) Vi = U^U^ 1 , 



then 



2 



(3.17) c = V 1 c : 

Repeating these steps results in a process given in (3.12) after an appropriate renum- 
bering of the c n . The n— th step of this process (|1.7p . results in the replacement 

(3.18) 



Vi -> 




C/i -> 




c/ 2 -> 


f/ 2 ,n 


£7 


On- 



4. Quasiperiodic System 

A crucial difficulty with perturbative schemes, is that these often produce terms 
with arbitrarily slow frequencies. Such terms lead, upon integration to small di- 
visors. Here, we will treat such a problem, and show, how to get the large time 
behavior of the system, to any order, via multi-scale time averaging. Our main 
example is the following: 

Example 1. Let Mj be N x N matrices with norm bounded by ||Mj|| < j^ 1 ^ 5 1 
for all j > 0. Assume, 

oo 

(4.1) A = Y^{Mjj W}i + c.c) 

where c.c stands for hermitian conjugate. Note that norm convergence of the sum 
is assured by our assumption on the norms of Mj. The case of interest is when 
ujj — > as j — > oo. 

In order to apply our method to this we have to split A (t) of (|4.1[) as 

(4.2) A = A>+A<, 

and to treat separately the two parts. For this purpose we introduce the following 
two Lemmas. 

Lemma 2. Let 

oo 

(4.3) A> = (Afj-e*"** + c.c) 

3=30 
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with jo sufficiently large, so that 

(4.4) UjT < 1, 

with T = /3~ 1//2 7 (3 small, and 



(4.5) 

for all j > j a . Then, 
(4.6) 

Proof. Start from , 
(4.7) 



\\M \\T Q <-^ 



f{A>{t')-A>\t'))dt' 
Jo 



S > 0, 



<0(1) 



and choose a term and denote it 

(4.8) B(t) = Mj-e^* 



(4.9) 
and 



(4.10) 

Define 
(4.11) 
for 

(4.12) 
Then, 

rt 



-mi 1 /" T ° e iw j T o - 1 

B o 0) = jr B(t)dt= e , 



J J n 



(n+l)T i(n+l)T Wj _ p inT u!j 
e*"'** = Af,- — — 

e iT Uj _ ^ s 



M,-e inT ° . 

iujjTo 



B s (t)=B { Q n) 
nT <t<(n + 1)T . 



dt' [B(t') - B 9 (t')] = / B(t')-B^(t') df + •••+/ B(t')-B^(t') 

JO Jo 1 J ijT L 



dt 7 



(4.13) +••• + 



where 

Since, 
(4.14) 
one finds, 
(4.15) 



'"To 

nTo < t < (n + 1)T . 

i r {i+1)Ta -r-i 

1 JjT 
rU+i)T 

/ [fl(t') - 5(0] ^' = 0. 

JjTa 
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dt'[B{t')-Bl(t')] = Mj 



iUjTg 



(4.16) 



71/ 



Je^t-nTo) - l) - (i-7iT ) r 



For, 

(4.17) Wj-To < 1, 

where T = ft- 1 ! 2 . Using the fact that \t - nT \ < T we get: 



(4.18) 

Assuming, 

(4.19) 



[B{t')-mt')]dt> 



\MA\T Q < 



< \\MjW -2T . 



J 



and summing the various contributions in the sum for A > , one finds, 



(4.20) 



[A>{t')-A>%t')]dt') 



< 



^2 2| jA^jjjTb < const < oo 



3=3a 



or 



(4.21) 



df (A>(i')-A> 9 (0) 



<0(1). 



Remark 2. Note the difference between ()4.2ip and f|2.50|) . 
Lemma 3. Let 

jo-l 



(4.22) 

and /ei 
(4.23) 

and T = Z?" 1 / 2 . T/ien 

(4.24) 

Proof. 

(4.25) 

we note that, 
(4.26) 



|M,-|| const 



'0 Jo 



A<(s)ds 



< const/3 



1/2 



30-1 



To 



□ 
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therefore, 
(4.27) 







< 



Therefore, the relevant sum is bounded as 

el 



(4.28) 

Therefore, 
(4.29) 



A> (*') dt' 



^ 2 II Af ,-|| 
< — — 3 — < const. 



OJj 



1 



A> (*') dt' 



< const/3 



1/2 



□ 



Remark 3. The partition (|4.2p is always possible since if (|4.4I) is not satisfied then 
ujjTo larger then some constant that is much smaller than unity. Therefore, all 
terms that do not belong to A > belong to ^4 < with the appropriate choice of the 
constant. 



Theorem 1. The global average A 9 2 defined by {2.5J$ is bounded by 
(4.30) |^|<0(^ 2 ) 



Proof. First we estimate A\. For this we define 
4 0) - ±f*°U^[A-Al]U Q 

U^{t)(f[A{s)-Al{s)] ds)u (t)] 



(4.31) 



1 

To 

i_ r T ° du-\t) 

T Jo dt 



A{s)-A Q 0) (s) dsU {t) ) dt 



To Jo 

From the definition of A °\s), 



A(s)-A^(s) 



ds I — -dt 

dt 



(4.32) 



ds 



(a(s)-A °\s) 



The contribution from the regime where Lemma [5] holds is, 



(4.33) 



ds [A(s)-A°(s)] 



<0(1) 



The contribution from the region where Lemma [3] is relevant is, 
(4.34) 



^— ^ t.tl 



e iUj(t-nTo) _ e in 



T _ ~ nT 0) 

To 



< 



□ 



\M* 
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and an equality similar to (|4.33[) holds. From the general theory (|2.10p and (|2.11ll . 

T~ 1 



(4.35) 



dU 



dt 



< O(p), 



dU n 



dt 



<0(f3) 



Therefore, combined with (|2.36l) 

(4.36) \Af\ < 0(j3) 

by (glU and (12~531 



.4 



(0) 



To 



1 

Tq Jo 
Jo 



1 



u[-\t') [AtW-AW)] U^dt' 

/•To 


" [A^s')- A\ (a')] ds') U!(s)ds 



[Ai(8)-A(\ d8\Ux{t) 



(4.37) 



ds 

r T 



k So U{1 (T [Al{s,) ~ A9i(s ' )] ds ') i Ui(s)ds - 



By (|2.50p of the general theory, 
(4.38) 



1 

To 



[A 1 (s)-Al( 8 )] dt 



<0(1). 



From the definition of A{ one finds f Q T ° [A x (s) - A{ (s)) ds = 0. Using (|2~37P 
combined with ([239]) and (|2~40| we find 



(4.39) 



d^ 



dt 



<0(f3 2 



dU~ 



dt 



<o(p 2 



then (|4~37l) combined with (|4~38l) and (f2T40|) leads to the bound, 



(4.40) 



A 



(0) 



< 0(/3 2 



But in each interval of length To is equal to one of the A% , each satisfying 
(|4.40l) . therefore it also satisfies (|4.40p . therefore 



(4.41) 



\A a A < 0([3 2 



Remark 4. This bound is better than the one of the general theory (12.541) . 
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